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ABSTRACT 

We report the discovery of a 4.5 /im counterpart to the anomalous X-ray pulsar (magnetar) 
IE 2259-1-586 with the Spitzer Space Telescope. The mid-infrared flux density is 6.3 ± 1.0 /xJy at 
4.5 /im and <20/iJy (at 95% confidence) at 8 /im, or 0.02% of the 2-10 keV X-ray flux (corrected for 
extinction). Combining our Spitzer measurements with previously published near- infrared data, we 
show that the overall infrared emission from IE 2259-1-586 is qualitatively similar to that from the 
magnetar 4U 0142-1-61. Therefore, the passive X-ray-heated dust disk model originally developed for 
4U 0142-1-61 might also apply to IE 2259-1-586. However, the IR data from this source can also be 
fitted by a simple power-law spectrum as might be expected from magnetospheric emission. 
Subject headings: infrared: stars — pulsars: individual: IE 2259-1-586 — stars: neutron — supernovae: 
general 



1. INTRODUCTION 

The notion of supernova fallback, where some of the 
ejecta from a core-collapse superno va ends up cap tured 
by the newly formed neutron star (|Chevalieij|1989[ ) and 
may have sufficient angular momen tum to form a disk 
()Lin. Wooslev. fc Bodenheimeilll99l[ ) , has been a general 
predi ction of supernova models (e.g.. I Wo oslev fc Weav er 
fl995h . Fallback can have profound effects on the final 
state of the ne utron star, forming disks that can give 
rise t o planets (jPhinnev fc Hansen! 119931 : iPodsiadlowskil 
Il993f ). and possibly even causing the neutron star to col- 
lapse into a black hole. Fa llback disks should manifest as 
a the rmal infrared excess (jPerna. Hernquist. fc NaravanI 
|2000D . but previous searches f or such excesses aro und 
neutron stars were unsuccessful (jLohmer et al.|[200l . and 
references therein). 

Several years ago, though, we discovered the 
mid-infrared (mid-IR; here, 4.5 and 8.0 /im) 

counterpart to the mag netar 4U 0142-f61 

(I Wang. Chakrabartv. fc KaplanI l2006l . hereafter 
IWCK06f l. The combined optical/IR spect rum of 
this magnetar suggests that the optical and IR data 
arise from two different spectral components. While 
the optical co mponent is demonstr ably of magneto- 
spheric origin (|Kern fc Martini l200l . we showed that 
the IR component ma y arise from a p assive (i.e., 
not accreting, but see lErtan et al.1 l2007f l. dust disk 
irradiated by X-rays from the ma gnetar (|WCK06l : 
I Wang. Chakrabartv. fc KaplanI l2008bD . The inferred 



spectral shape, while not well constrained, is remarkably 
simila r to those of protopl anetary disks around young 
stars (i Beckwith et al.iri990D . The disk's survival lifetime 
10^ yr) significantly exceeds the pulsar's spin-down 
S'gs'^ (^ 10^ yr), consistent with a supernova fallback 
origin. 

This dust disk model, while intriguing, is not definitive. 
In particular, we have yet to firmly establish whether the 
mid-IR flux in 4U 0142-f 61 comes from a disk or from the 
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gle wavelength or across wavelengths can be a powerful 
discriminant , depending on t he t imescale. This was il- 
lustrated bv lKern fc Martini (|2002i) . who found that the 
pulsed fraction in the optical exceeded that in soft X- 
rays, making it impossible for t he optical to res ult from 
reprocessed X-rays (also see Dhillon et al.ll2005[ ). In the 
infrared the si tuation is less clear. Th e significant vari- 
ability seen bv lDurant fc van KerkwijS (j2006Q) at 2.1 ^m 
with no related change in soft X-rays pro vides a stiff test 
for a disk model, but is not definitive fsee lWang fc Kaspl 
l2008f ). The unknown behavior at higher energies, the 
complicated behavior across the optical/near- infrared 
(near-IR), and the lack of pulsations at 2.1 /xm (< 17% 
pulsed at 90% confidence: iMorii et al.1 [20091 ) all compli- 
cate the matter. 

As we attempt to determine the origin of the mid- 
IR emission, we must equally attempt to understand 
how ubiquitous mid-IR counterparts to young neu- 

^ While spin-down is not an accurate measure for the age of 
a magnetar ( Woods & Thompson 2006), no independent age esti- 
mate is available for 4U 0142+61 and this age is consistent (if not 
somewhat longer than) typical ages of other magnetars. 
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tron stars are. Suggestively, all four of the mag- 
netars with confirmed quiescent near-IR counterparts 
have the same near-IR/X-ray flux ratio of « 10~^ 
(jPurant fc van KerkwiiMl2005ll2006"d[) . In the disk sce- 
nario, this ratio should be determined largely by geom- 
etry, so mid-IR observations of other magnetars down 
to this fractional level may help establish whether disks 
are present, but similar behavior might also be expected 
from magnetospheric emission. Shallow mid-lR searches 
of three other magneta rs we re not very constraining 
(|Wang. Kaspi. fc H i gdon|[2007t ). but deep a 4.5 fj,m upper 
limit for IE 1048.1—5937 following an X-ray flare con- 
cluded that it did not have a mid-IR counterpart similar 
to that of 4U 0142-1-61 despite its similar ncar-IR/X-ray 
flux ratio, thus casting some doubt on the disk interpre- 
tation (Wang et al. 2008a). Here, we present the results 
of a deep search with the Spitzer Space Telescope for a 
mid-IR counterpart to the magnetar IE 2259-1-586. 

IE 2259-1-586 was identified as an X-ray point source 
with coherent 7 s pulsations in the center of the su- 
perno va remnant (SNR) CTB 109 ( Fahlman fc Gregorvl 
|1981|) . The association leads to a distance of 3.0±0.5 kpc 
based on interactions betwee n the SNR and H II re- 
gions with meas ured distances (iKothes. Uyaniker fc Yail 
|2002|) . although IFurant fc van KerkwiikI (|2006a[ ) deter- 
mined a distance to IE 2259-1- 586 of 7.5 ± 1.0 kpc us- 
ing the "red clump" method. iDurant fc van KerkwiikI 
()2006bD found the column density to IE 2259-^586 to 
be TVh = (1.1 ± 0.3) X 10^^ cm^^, corresponding to a 
visual extinction of Ay = 6.3 ± 1.8 mag, from fitting 
of X-ray a bsorption edges (' consistent with earlier mea- 
surements; iPatel et al.ll200l . The near-IR (ii's-band, or 
2.1 jLtm) counterpart t o IE 2259-1-586 was identified by 
iHuUeman et al.l (|200lD . Uniquely among magnetars, the 
near-IR flux was observed to vary in concert with the 
X-ray flux, with both declining following a major series 
of X-ray bursts fFam et al.ll2004l ). This led to the sug- 
gestio n of X-ray heating of a disk (as in iMatsuoka et al.l 
Il98l for a low -mass X-ray binary), making this source 
an especially attractive target for mid-IR searches. 

2. OBSERVATIONS AND REDUCTION 

2.1. IRAC Data 

Our primary data were observations of I E 2259-1-586 
with the Infrared Array Camera (IRAC; 'Fazi o et al.l 
2004 ) onboard the Spitzer Space Telescope ( Werner et al.l 
2004li on 2007 August 8. We observed IE 2259-^586 us- 
ing two of the four possible IRAC channels: 4.5 ^m and 
8.0 /im (IRAC channels 2 and 4). The observation con- 
sisted of 50 dithered exposures of 96.8 s for a total inte- 
gration of 80.7 minutes. 

We started with the Basic Calibrated Data (BCD, from 
pipeline version S16.1.0), discarding the first exposure at 
each wavelength as recommended. We then processed 
the BCD images though an artifact mitigation pipeline^ 
(2005 October 13 version). With the MOPEX (MOsaicker 
and Point source Extractor, ver. 16.3.7) pipeline, we mo- 
saiced the individual exposures together by interpolating 
the data onto a common grid and rejecting radiation hits. 
The final images are shown in Figure [T] 

2.2. Gemini Data 



Our analysis requires precise relative astrometry, and 
the angular resolution of the IRAC images results in 
many blended objects. To aid in interpreting our mid- 
IR data and in particular to serve as an improved as- 
trometric reference, we analyzed an archival near-IR 
{Ks-hand, or 2.1 fim) o bservation taken wit h the Near- 
Infrared Imager (NIRI; iHodapp et aIll2003D on the 8 m 
Gemini North telescope at Mauna Kea, Hawaii. The 
data consist of 50 exposures with the //6 ca mera taken 
on 20 03 May 27 (and previously published bv lTam et al.l 
|2004|) . each with 4 x 15 s integrations, for a total expo- 
sure of 50 minutes. We used an IRAF package available 
from the NIRI Web site^ to reduce the data. We pro- 
ceeded through the steps of this package, flatfielding the 
data, subtracting the sky, shifting the images, and adding 
them together. We referenced the astrometry of the fi- 
nal image to the Tw o Micron All Sky Survey (2MASS; 
ISkrutskie et~alll2006[ ). finding 57 stars that were not sat- 
urated or badly blended, and getting rms residuals of 
0'.'16 in each coordinate. This image is also shown in 
Figure [TJ 

3. ANALYSIS 

We searched the IRAC images in Figure [T] for a 
counterpart to IE 2259-1-586 by looking for a point 
source at the position correspo nding to star 1 (which is 
IE 2259+586; HullejnaneLiD [2001) in the NIRI image. 
At this position, we see a faint point source in the IRAC 
4.5 /Ltm image (hereafter the "IRAC source"). However, 
one must be careful, as the position of star 2 only dif- 
fers from that of the magnetar by 0'.'97; this is «1 IRAC 
pixel width, and only half of the angular resolution of 
the image (« 2" FWHM). 

We performed photometry on the 4.5 /im IRAC im- 
age to find the positions and flux densities of all of the 
point sources. We used the APEX (Astronomical Point 
Source Extraction, part of MDPEX) software to iden- 
tify and perform point-response function (PRF) fitting 
for the photometry using the routine apex_lfrcLme.pl 
with the Spitzer-snpplied PRF. The PRF is the tradi- 
tional point-spread function (PSF) convolved with the 
pixel-response function to determine how a source ac- 
tually appears in the data^. Part of this analysis iden- 
tified blended objects using simultaneous PRF fitting, 
but the IRAC counterpart was consistent with a single 
point source (the per degree of freedom for the PRF 
fit was 0.48, consistent with other point sources in the 
field which have reduced of 0.1-3) with flux density 
6.3±1.0 /iJy at 4.5 ^m. The uncertainty here includes the 
standard term from the PRF fitting but is dominated by 
substantial contributions from difficulty in robustly mea- 
suring the background and identifying which pixels make 
up the source. At 8.0 /im there is no detection, and we 
estimate a limit of < 20 /iJy (95% confidence, limited 
largely by the varying diffuse background) based on the 
APEX detections in that region. 

In determining whether the IRAC source is indeed the 
counterpart to star 1, we considered three issues. First, 
we were concerned with whether the position of the IRAC 
source is consistent with that of star 1 or 2. We took 

^ See http : //www.us-gemini .noao . edu/sciops/instruments/niri/NIRIIndex 

^ The distinction between the PSF and PRF is important 
for undersampled or critically sa mpled dat a like those here; see 



See |http: //spider . Ipac ■ caltech. edu/staf f /carey/irac_artif act^b^tp : //ssc . spitzer . caltech.edu / irac/ps f .html | 
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Fig. 1.— Field of IE 22-59+586 at 2.1 /im (NIRI, top left), 4.5 ^tm (IRAC, bottom left), and 8.0 /^m (IRAC, top right). On the IRAC 
images, we label several astrometric reference stars (A, D, and F) from iHuUeman et al.l ([2001') and from this work (101 and 102), along 
with the IRAC source. In the NIRI image, we label the same reference sources, plus we label stars 1 and 2 separately. In the bottom right, 
we show a zoom around the position of IE 2259+586 in the 4.5 /im IRAC image (the region is indicated by the box in the larger image). 
The position of the IRAC source is shown by the circle, with the circle's radius equal to the astrometric uncertainty (0'.'15). The positions 
of stars 1 and 2 from the NIRI image are also shown with appropriate error circles. In all images, north is up, east is to the left, and scale 
bars are indicated in the lower left. 

the astrometry (using both APEX and standard centroid- 
ing; both gave consistent resuhs) of the eight reference 
sources (some of which are labeled in Figure [1]) , along 
with the NIRI astrometry of those objects and stars 1 
and 2. We used the reference stars to refine the astrom- 
etry of the IRAC image relative to the NIRI image (we 
did not use a proper reference frame as an intermediary, 
as this would have introduced additional uncertainties). 
Fitting only for an offset (the position angles of both ob- 
servations were referenced independently to 2MASS), we 
found a small shift (O'.'l) with an rms of 0'.'05. With this 
shift, the position of the IRAC source differs from that of 
star 1 by 0'.'28, and from that of star 2 by 0'.'74 (see Fig- 
ure [2|) . There is an intrinsic centroiding uncertainty for 
the IRAC source (the potential counterpart) of k, 0'.'15 in 
each coordinate (90% confidence radius of 0'.'32). Given 
this, the IRAC source is largely consistent with the posi- 



tion of star 1 (the probability of chance alignment is 78%) 
and is inconsistent with the position of star 2 (the prob- 
ability of chance alignment is 99.998%). The situation 
is also illustrated in the bottom right panel of Figure [1] 
which shows the extent of the PRF of the IRAC coun- 
terpart, the measured centroid, as well as the positions 
of stars 1 and 2 in the same reference frame. 

Second, we compared the IR colors of the IRAC source 
to all of the other objects in the field. We show a color- 
color diagram in Figure [3l For the IRAC source, we 
used the near-IR (J- and K^^h&nd) photometry of star 1 
given by iHulleman et all (|2001D: for star 2 we also used 
the photometrv from lHulleman et al.l (|2001h ; for the field 
stars, we used 2MASS for the near-IR. Assuming, as de- 
duced from the above astrometry, that the IRAC source 
is associated with star 1 and not with star 2, we see that 
star 2 is consistent with the bulk of the field population, 
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Fig. 2. — Relative astrometry of the IRAC source. We show the 
residual position difference for eight reference sources (black points) 
that we used to transform between the IRAC image and the much 
higher resolution NIRI image after accounting for a net shift of 
O'.'l. We also show the measured position of the IRAC source (blue 
circle) and the positions of star 1 (red diamond, labeled "AXP" ) 
and star 2 (red square). The position of the IRAC source is much 
more consistent with that of star 1 than star 2. Also see the bottom 
right panel of Figure [T] 

i.e., reddened main-sequence stars. On the other hand, 
the IRAC source is significantly redder in both J — Kg 
and Kg — [4.5 /Ltm] and clearly stands out from the popu- 
lation. It is quite close in colors to 4U 0142-1-61, for the 
same extinction. Uncertainties in the visual extinctions 
of both IE 2259-^58 6 and 4U 0142-1-61 amount to ±2 mag 
(jPurant fc van Kerk wiik 2 GG6b) and do not change the 
general agreement because of the small extinctions at 
near- and mid-IR bands {Aj « O.SAy, Ak ~ O.l^y, 
^4.5/im ~ O.OSAv^; also see Figure S])). 

Finally, we examined the spectral energy distributions 
(SEDs) of the IRAC source and star 2 (Figure H]). First, 
we estimated the spectral type of star 2. With only 
two data points (J = 23.1 ± 0.1, = 21.5 ± 0.2) we 
could not discriminate between later stellar type and 
higher reddening (the H- and /-b and upper limits from 
iHuUeman. van Kerkwiik. fc Kulk arni 2000 were not con- 
straining), but we find reasonable solutions for K/M 
dwarf stars with a few magnitu des of extin ction and dis- 
tances of a few kpc (based on ICox|[200Clf ). For stars of 
those spectral types, we would expect the SED to peak 
in the near-IR region and decline in the IRAC bands, 
with typical flux densities of ~ 0.5 /zJy expected. Even 
brown dwarfs, which can be redder in JtTg — [4 .5 /.tm], have 
F^(2.1 ^m) > i^^(4.5 Atm) (jPatten et al.ll2006t l. There are 
classes of objects such as young stellar objects or massive 
stars with dusty winds that could produce suc h red colors 
(e.g.. lUzpen et al.ll2007t iMuench et"ani2007f) . but they 
are relatively rare (although this is a comp licated region 
with many objects along the hue of sight; iKothes et al.l 
l2002f ). often brighter than star 2 (spectral types earlier 
than about FGV would have to be well outside the Galaxy 
to have the observed Ks magnitude), and would likely 
have been identified at this position in other wavelengths 
(X- ray, near-IR, molecular line, or continu um radio, etc.; 
see IKothes et al.l 120021 : ISasaki et al.ll2004[ ) but were not. 
In contrast, the 4.5 \xm flux density of the IRAC source 
is a factor of 3 higher than the /sTg-band flux densities of 
star 1 or 2, and an order of magnitude higher than the 
expectation for a normal star at 4.5 ^m. 

Taken together with the astrometry, this is strong sup- 
port for not associating the IRAC source with star 2, 
but instead with star 1 . It is extremely unlikely that the 
IRAC source is associated with neither object, as proba- 
bility of chance coincidence with star 1 is only 0.3%, and 
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Fig. 3. — J — color vs. — [4.5/im] color for the field of 
IE 2259+586. The objects from the field with both 2MASS and 
IRAC detections are the black points. IE 2259+586 is the upper 
limit in the upper right, while star 2 is also an upper limit (this up- 
per limit is only approximate, since star 2 is located within the PRE 
from IE 2 259+586). We also plot the co unterpart of 4U 0142+61 
IIWCK06I ) reddened from A^/ = 3.5 mag HDurant fc van KerkwiiH 
I2006bl) to the same total extinction as IE 2259+586. A redden- 
ing vector for Ay = 5 is shown, as is a main seq uence with zero 
reddening (solid line, computed from IKurucd[T993l models: a giant 
branch will appear very similar). 

this would also imply even more extreme colors (Fig- 
ure [3|) that would be clearly non-stellar. Even with sev- 
eral tens of magnitudes of extinction it would be dif- 
ficult to move the IRAC source onto the stellar locus, 
and examining the NIRI image we see no other sources 
with j^ s — [4.5 /im] > 2. From [Purant fc van KerkwiikI 
(|2006al) . an extinction of only ~ 10 mag is reached at 
~ 8 kpc, and the molecular gas maps of Colden^ find a 
maximum extinction of ^ 5 mag, so having very high ex- 
tinction outside of very local extinctions seems unlikely. 
We therefore conclude that the IRAC source is indeed 
the mid-IR counterpart to star 1/lE 2259+586. 

4. DISCUSSION 

We have found a 4.5 /im counterpart to the magne- 
tar IE 2259-1-586, with flux density 6.3 ± 1.0 ^Jy and a 
limit of < 20/iJy at 8.0 /im. The absorption-corrected 
mid-IR/X-ray® flux ratio for this source is 1.9 x 10"''^, 
very similar to the ratio for 4 U 0142-1-61, so the com- 
mon X-ray to 2.1 /im flux ratios (jPurant fc van Kerkwiifl 
[2005h seem to extend further into the infrared. We 
note, though, that one must be careful with any model 
or interpretation for the IR emission since the 2.1 
flux density of IE 2259+586 (and of 4U 0142-^61 for 
that matter ; iHuUeman, van Ker kwiik. fc Kulkarni|[200l 
iDurant ^F^ an Kcrkwijk 2006c) is known to vary. The 
2.1 fim flux density we plot in Figure |4] is the faintest 
measured and is a ssumed to be close to the baseline level 
(|Tam et al.ll2004D and should not be affected by flares, 
but it was measured two years before the IRAC observa- 
tions. For IE 2259-1-586, the I R variability had been 
seen along with X-ray flaring (K aspi fc Gavriill l2002t 
iKaspi et al.ll2003i) in the past. Since then, though, no 
transient X-ray behavior for IE 2259-f 586 was reported 
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See |http : //asc. harvard. edu/toolklt/colden. jsp[ 

8 The unabsorbe d X-ray flux of IE 2259+586 
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' ergs~^ cm~^ (Patel et alji 200 T). although this is only over 
the 2— lOkeV range and does not include softer emission which is 
uncertain due to the moderate absorption. 
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Fig. 4.— SED of IE 2259+586 and star 2. We plot uF^ (left axis) 
vs. frequency for IE 2259+586 (solid red points/upper limits) and 
star 2 (blue squares), and also give the luminosity vLi, (right axis) 
for a distance of Sda kpc. For star 2 we also show one of the possible 
stellar fits (an M5V star at Ay = 3). The open red points/upper 
limits are the data on I E 2259+586 deredd ened with Ay = 6.3 us- 
ing the reddening law of llndebetou w et al.' (2005), and the shaded 
band shows the ±1 cr range of Ay (Durant & van Kerkwiik 2006b). 
We plot a power law, Fi, oc , which fits both detections for 

IE 2259+586 (dashed line) and an irradiated disk model (based on 
IWCK0& : dot-dashed line). The various bands are labeled. 

and the X-ray flux is now near the quiescen t level (al- 
thoug h with a slightly different spectrum; IZhu et al] 
120081) . so our assumption of a baseline flux seems safe, 
but X-ray monitoring is not very regular and there could 
potentially be IR variability independent of X-ray activ- 
ity. 

With the limited data we have, a wide range of model 
fits are possible. Here we discuss two general categories: 
a power law and an irradiated disk model. A simple 
power law {Fi, oc v^^-^^) can fit the two detected IR 
flux points and is also below the upper limits at 1.2 /im 
and 8 /im (Figure |4]). Such a power law spectrum could 
arise fro m the pulsar magn etosphere (as for Crab pul- 
sar; see iTemim et al.l I2006D . which could also possibly 
produce correlated X-ray and infrared flux changes, al- 
though this has not been studied in detail. However, 
unlike the Crab pulsar, this power-law rises further into 
the mid-IR band, and it would be interesting to obtain 
longer- wavelength data on IE 2259+586 to see if the SED 
keeps rising past 4.5 /im, but source confusion makes 
this very difficult. The rising spectrum is problematic 
for some detailed ma gne tospheric mod els (like those of 
lErtan fc Chen3i200l and iHevl fc HernquistH 2005l . but a 
generic magnetospheric origin is c ertainly p ossible. 

However, as with 4U 0142+61 ()WCK06l ). we can also 
fit the IE 2259+586 d ata with an irradiated passive disk 
model (also see Vrtil ek et all [1990; Pcrna ct al. 2000; 
other spectral shapes are also possible, but the fits are too 
unconstrained and without specific motivation we will 
not address them) . Without more data, we cannot make 
conclusions about the presence or absence of a disk, and 
indeed other observa tions suggest that a disk interpreta- 
tion is problematic (|Wang et al.l l2008al) . However, it is 
still possible, and we fi nd it worthwhile to expand on the 
discussion of IWCK06I and further explore the implica- 
tions such a disk would have, noting in particular where 
IE 2259+586 differs from 4U 0142+61. 

This model fit is qualitatively similar to what we ob- 
tained in WCK06 for 4U 0142+61, with a smaU inner 
radius and an outer radius a factor of a few to 10 larger. 



The details of both fits depend on assumptions about the 
distances, extinctions, and inclinations, but the rough 
shapes of the SEDs going from the near- to mid-IR 
are similar, with a factor of 3 in fiux density increase 
from 2.1 to 4.5 /xm, after co rrecting for extinction (un- 
like what IWang eral]l2008al found for IE 1048.1-5937). 
The inner radius is constrained by the need to have 
the -fCg-band measurement lie below the 4.5 /im mea- 
surement, while the outer radius is less constrained by 
the 8.0 /im upper limit. Note that the nominal inner 
radius of i?in = 0.25i?o (for a distance of 3 kpc and 
an inclination of 60°) is smaller than the light cylinder 
radius (O.4i?0) which would imply int eraction between 
the d isk and the magnetar's spin (e.g.. IChatteriee et al.l 
12000) — possibly contradicted by its observed steady spin- 
down (Gavriil fc Kasp i 2002 ). This could either point to 
an inconsistency in our model, or simply that we must 
choose the other model parameters (distance, inclina- 
tion, etc.) such that the inner radius is > O.5i?0; as 
an example, a larger dis tance such as that determined 
bv lDurant fc van Kerkwi ik (2006i|) has an inner radius^ 
more than twice the light cylinder radius. 

The b inding energy of the putative disk, following 
IWCK06I to estimate the disk mass, is ~ 10''^ erg for the 
upper end of the mass estimate. This is large enough that 
we do not have to worry that such a disk would have been 
unbound by the large X-ray bursts (a total fluence of 
3 X 10~^ e rg s~^ , or total isotr opic energy of 3 x 10^ ''^3 erg) 
found bv lKaspi et al.l ()2003f ). and the physical state and 
structure of the disk may not be altered dramatically, 
as this energy is equivalent to only ~ 10'^ s of normal 
X-ray activity. Even a disk several order of magnitude 
lower in mass would remain bound. Only a truly giant 
flare, such as those seen from several soft 7-ray repeaters 
(| Woods fc Thompson|[2006l ). would be energetic enough 
to disrupt a disk like this, althou gh the details o f that 
process are difficult (also see iWachter et al.ll2008| ). but 
we stress that our mass estimate is an upper limit, and 
at the low end of the allowed mass range the disk would 
not be bound. Again, this could indicate either a problem 
with the disk model, or that objects with disks cannot 
have had large flares. 

Among magnetar s, the corre l ated n ear-IR and X-ray 
variability found by iTam et al.l ()2004[ ) for IE 2259+586 
is unique. Such a correlation is a natural (but not neces- 
sarily unique) consequence of the disk model: a change in 
the X-ray flux would produce an accompanying change 
in the infrared flux from the disk, and while the to- 
tal reprocessed flux is a constant fraction of the X-ray 
flux, the IR flux at a given wavelength can change faster 
or slower as an nuli of a given temperature mov e and 
change area (cf. Ivan Paradiis fc McClintockllT99l. We 
find t hat the 2.1 /im flux increases seen by iTam et al.l 
(l200l . which were of comparable amplitude to the X- 
ray flux increases, are roughly consistent with a disk like 
that considered here, although this assumes that 100% 
of the 2.1 /im flux comes from the disk and ignores the 
possibility of physical changes to the disk such as move- 
ment of the inner radius following the flux changes. How- 

^ As mentioned in IWCK0(3 . the inner radius discussed here is 
just that to which dust can penetrate. Gas could extend further 
inward, so that even with the inner dust radius larger than the 
light cylinder radius interactions could still occur. 
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ever, correlated changes were n ot seen in 4U 0142+61 
(jPurant fc van Kerkwiild[2006d ). who instead saw rela- 
tively rapid, significant changes at 2.1 /xm without any 
changes in X-rays. Any model must be able to accom- 
moda te this wide range in behavior (see IWang fc Kaspil 
|2008| ). Additional observations at longer wavelengths 
where the decomposition is less ambiguous could again 
provide important constraints. 
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